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Critical Exponents of the Classical Heisen- 
berg Ferromagnet 

In a recent letter, Brown and Ciftan (BC) [1] reported 
high precision Monte Carlo (MC) estimates of the static 
critical exponents of the classical 3D Heisenberg model. 
While their finite-size scaling (FSS) analysis yields values 
for the critical temperature Tc = l/Kc and the critical 
exponent ratios l3/i>,'j/i>, which are compatible with all 
recent findings, BC claim that the specific heat C of this 
model is divergent at Tc, which is in strong disagree- 
ment with other recent high statistics MC simulations, 
high-temperature series analyses, field theoretic methods 
[2-5], and experimental studies [6], which all find a finite 
cusp-like behavior. 

In their Ansatz (2) BC use a non-linear six parameter 
fit to 14 data points for C on lattices of linear size 2 < 
L < 32. The fit resulting in a/iy = 0.117(4) has still 
a total ^ 60, and therefore is by standard reasoning 
not acceptable. This indicates that either the statistical 
errors of their data are underestimated or the use of such 
small lattices as _L = 2 requires the inclusion of even 
more correction terms. Alternatively one may ask how 
the fit parameters would change if the smallest lattices 
are successively discarded. By hyperscaling BC deduce 
from this value a non-standard exponent i> = 0.642(2), 
leading in turn to non-standard estimates of /3 and 7. 

We find it very dangerous to base such an incisive con- 
clusion solely on the very delicate FSS behavior of the 
specific heat. In particular we strongly disagree with the 
statement of BC that i> is extremely difficult to mea- 
sure directly. The derivative of the Binder parameter 
dU I dK , and the logarithmic derivatives d\n{m) / dK and 
dlnirn^) I dK all scale like L^l'^ and, using fiuctuation for- 
mulas, can be as easily measured as C , and the statistical 
errors are straightforward to control. Already our data 
for dU/dK in Ref. [2] gave an estimate of = 0.704(6) 
which agrees with our estimate of a/i> by hyperscaling. 
Moreover it is compatible with the value i> = 0.698(2) de- 
rived from fits to the critical behavior of independent cor- 
relation length data in the high-temperature phase [2]. In 
Ref. [4] we studied this model with emphasis on topolog- 
ical excitations on much larger lattices with 8 < L < 80. 
By analyzing the new data with the above three quan- 
tities at K = 0.6930 Kc we obtain from FSS fits a 
prediction of = 0.699(3) (cp. Fig. 1). 

For our large lattices of Ref. [4] a fit of the cusp 
form C{L) = C'^s + CqL"^" to 12 data points yielded 
a/iy = —0.225(80) (x^ = 7.84). A more precise estimate 
can be obtained by analyzing the energy directly accord- 
ing to E{L) = E'^'^s + £'o_L(«-i)/'', which definitely has a 
regular term. This yielded again a negative exponent ra- 
tio a/v = —0.166(31) (x^ = 11.3), which by hyperscaling 
is completely consistent with our value of i^. If we allow 
for a confiuent correction EiL'^"^~^^^'^ , then the fit only 
slightly improves (x^ = 10.9) with an almost vanishing 
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FIG. 1. Double logarithmic plots of data and fits at 
K = 0.6930 for dU/dK (yielding = 0.703(6)), d\n(m)/dK 
(ly = 0.698(3)), and dln{m'^}/dK (i^ = 0.697(3)) versus L. 

amplitude Ei = — 0.8x 10""'^^, showing that there really is 
no need for an additional term. If we add a term CiL"^^'^ 
to the C-fit, again we only get a marginally improved fit 
(x^ = 7.79), but now with a positive exponent, a/v = 
0.09(1), and ai/i> = —0.6(1). However, there are also 
different solutions with about the same x^, showing the 
danger of being misled by a too fiexible C-fit Ansatz. 

The "universality scaling" of Fig. 3 in [1] indeed looks 
best for i> = 0.642, but this plot cannot serve as an in- 
dependent determination of i>. The predicted data col- 
lapse is an asymptotic statement for large L near Tc, and 
neither should it be expected far beyond Tc, nor for very 
small lattice sizes, even if the correct exponents are used. 

The numerical coincidence that Kc ^ lii(2) should be 
interpreted with extreme care and by no means implies 
that the model is exactly solvable. Even for the simpler 
3D Ising model no analytical solution is known - and the 
6th digit of a similar conjecture for Kc has recently been 
questioned numerically [7]. 

To summarize, we see no compelling reason that one 
should have any doubts that the estimates of the critical 
exponents of the 3D Heisenberg model, as consistently 
obtained by field theory, high-temperature series expan- 
sions, and MC methods, need to be reconsidered. 

Christian Holm 

Institut fiir Theoretische Physik, FU Berlin 
Arnimallee 14, 14195 Berlin, Germany 

Wolfhard Janke 

Institut fiir Physik, Johannes Gutenberg-Universitat 
Staudinger Weg 7, 55099 Mainz, Germany 


April 19, 1996 
PACS numbers: 


75.10.Jm, 05.30.-d 


[1] R.G.Brown and M. Ciftan, Phys. Rev. Lett. 76, 1352 (1996). 

[2] C. Holm and W. Janke, Phys. Rev. B48, 936 (1993). 

[3] K. Chen et at, Phys. Rev. B48, 3249 (1993). 

[4] C. Holm and W. Janke, J. Phys. A27, 2553 (1994). 

[5] J. Adler et at, Physica A201, 581 (1993). 

[6] S.N. Kaul, J. Magn. Magn. Mater. 53, 5 (1985). 

[7] H.W.J. B15te et al, J. Phys. A28, 6289 (1995). 


1 


